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ABSTRACT 


Interest in radar imaging has been growing for the past decades because of its 
long range sensing capabilities. It’s continued utility and applications in wide-ranging 
areas is fundamentally dependant on the ability to produce high-quality, artifact-free 
imagery. The use of radar to identify and image moving targets remains of great interest 
for both commercial and military usage. However, when imaging moving targets, there 
will be issues of incorrect positioning or streaking, as the unknown target velocity gives 
rise to image artifacts. Many techniques have been developed to handle moving objects, 
however, these techniques make use of the so-called start-stop approximation in which 
target motion is assumed to be momentarily stationary, while it is being interrogated by a 
radar pulse. 

A new linearized imaging theory that combines spatial, temporal and spectral 
aspects of scattered waves has been developed. This thesis considers the performance of 
these techniques and compared to existing imaging schemes. It is shown that the new 
imaging scheme provides better localization, and is translation invariant in phase-space. 
It is also shown that the imaging scheme is dependent on the aperture geometry. 
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I. 


INTRODUCTION 


A. RADAR IMAGING HISTORY 

Radar imaging is a technology that has been developed mainly within the 
engineering community. Radar is a sensor system that is widely used to detect, locate and 
identify targets at great distances in all kinds of weather. Radar’s ability to perform such 
functions unlike many optical sensor systems has made it a proven and well-recognized 
system for the past 70 years. The fundamental reason for radar’s usefulness is its ability 
to interpret and extract information from the echo signals, and the fact that wavelengths 
of radar signals makes them relatively unaffected by atmospheric and weather-induced 
attenuation. However, this behavior causes problems when using radar for imaging: 
Image resolution depends on signal wavelength. Radar systems generally use longer 
wavelengths and hence suffer reduced image resolution compared with other systems. 
Nevertheless, radar is still sensitive to objects whose length scales range from centimeters 
to meters. This is because radar waves scatter from objects whose size is on the same 
order as the transmitted signal wavelength or larger, and many objects of interest are in 
this range. 

Radar programs started in the 1930s and much research and development work 
has been done to improve the use of radar for imaging. Many radar imaging techniques 
have also been developed to increase the resolution of radar-based imaging. Examples of 
the key inventions are the Synthetic-Aperture Radar (SAR) and Inverse Synthetic- 
Aperture Radar (ISAR), which were developed due to interest in operational concepts. 
The key goal in both of these applications is the ability to produce high-quality, artifact- 
free imagery. 

B. IMAGING TECHNIQUES TO CREATE ARTIFACT FREE IMAGERY 

Radar-based imaging techniques have evolved considerably over many decades. 
One common objective among these developments is to obtain an artifact free image. 
Many implicit assumptions are embedded in these techniques which are discussed in 
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many journal papers. Among these assumptions, the start-stop approximation is common 
to all existing imaging techniques and this is central to the thesis and holds the 
fundamental difference between the new imaging scheme developed in [2] and current 
techniques. 

1. Linearized Scattering Model by Means of the Born Approximation 

Scattering models in the past have been built upon a foundation that views targets 
as being composed of simple, non-interacting point scatters. When the imaging scene 
consists of multiple targets, multiple scattering occurs and accounting for this makes the 
problem extremely complex. For this reason, the “weak scatterer” or Born approximation 
is used. This approximation is very useful because it makes the imaging problem linear 
and manageable. 

However, the Born approximation has its limitations that may result in image 
artifacts due to its neglect of the multiple-scattering terms. Actual multiple-scattered 
waves return to the radar later than expected from the single-scattering model. These 
waves are then interpreted as having come from an object farther from the radar. In 
addition, increased capabilities in radar resolution and sensitivity and parallel 
advancements in computational technology have exposed the inadequacies of such simple 
scattering models. With the ability of radar systems to resolve and image the components 
(subscatterers) that lie within the support of a traditional target, it is inevitable that new 
radar target models that can describe these multiple scattering effects accurately, whether 
multiple targets or multiple components of a single target are required. Current 
approaches typically utilize correction techniques to enhance imaging algorithms based 
on the low order “weak scatterer” model. 

2. Bandwidth Limited Radar Systems 

All real systems are bandwidth limited. In addition, radar data is often corrupted 
by noise. Consequently these limitations induce unwanted image artifacts and reduce 
image resolution. 
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3. 


“Start-Stop” Approximation 


In practice, modern pulsed radar systems do not usually measure the Doppler 
shift. Instead they estimate velocity from the phase difference between returns from 
successive pulses. Hence, the waveform commonly used by modem radars is a train of 
high-range resolution pulses that enable pulse-to-pulse velocity estimation. Such 
waveforms inherently invoke a start-stop approximation, which assumes that the target is 
stationary during the time the radar pulse illuminates the target and moves only between 
pulses. Because the individual pulses are short and targets of interest are small and move 
slowly compared to the speed of light, this is almost always a good approximation. Using 
this approximation, good velocity estimates are based on range-rate (two accurate range 
measurements, separated by a time interval). However, no Doppler shift measurements 
are actually measured by the system. From the measured range and range-rate data, 
spatial images are then formed. 

However, in practice there is a limit to the how short a pulse can go. In order to 
achieve a relatively good signal-to-noise (SNR) ratio by coherent integration, long pulse 
trains are required. As a result, when the returned pulses are coherently processed at the 
receiver, the stationary target assumption may no longer hold. The scatterer may move to 
a different position within that pulse interval or resolution cell and result in a blurred 
image. 


C. MOTIVATION-IMAGING MOVING TARGETS 

To date, there have been a number of attempts to develop imaging techniques that 
can handle moving targets. For example space-time adaptive processing (STAP) is a 
signal processing technique that was originally developed for detecting slow-moving 
targets via airborne radars. Multiple-element antenna array together with real-aperture 
imaging techniques to produce ground moving target indicator (GMTI) images [24] were 
the key features in STAP. The technique described in [25] together with SAR (designed 
to image stationary scenes) and GMTI processing develop a new techniques for detecting 
slow-moving surface targets that exhibit start-stop maneuvers. Another method described 

in [12] discusses the extension to non-sideways looking array radar. 
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The concept of velocity synthetic aperture radar (VSAR) is another development 
described in [10]. VSAR is essentially a multi-element SAR system involving 
conventional processing to form an image at each element. Here, the image phases are 
preserved and compared to deduce target velocity. However, VSAR processing assumes 
that the scatterers remain in a single resolution cell throughout the integration time. 

Time-frequency methods are also very useful in the area of image generation. 
This method, when used for signal analysis, allows decomposition of the frequency 
spectrum of time-varying signals in shorter time frames and provides an added dimension 
for examining the dynamic behavior of the signal as it varies over time. The size and 
shape of the weighting window function can be altered to fit the specific analysis 
requirement for the signal. Recently, this method was combined with ISAR to improve 
imaging, motion compensation and micro-Doppler target vibration studies [8 20 and 21]. 
It was reported in [20] that the time-frequency transform (TFT) was used to construct 
inverse synthetic aperture radar (ISAR) images of targets with high rotation rate by 
helping to overcome the migration of individual scattering points from one range cell to 
another. 

However, all these imaging techniques make a common approximation that a 
target in motion is assumed to be momentarily stationary while interrogated by a radar 
pulse. Therefore, it is apparent that the need for an imaging scheme that can 
accommodate target motion during imaging process is necessary and this constitutes a 
main focus of this thesis. 

D. OBJECTIVE 

The objective of this project is to study a linearized imaging theory developed by 
Professors Cheney and Borden [2] for an imaging scene containing moving targets and 
this will be done in four parts. Part 1 analyses the Point Spread Function (PSF) focusing 
on imaging a static point scatterer. This was completed in a previous thesis [3]. Part 2, 
which is the main focus of this thesis, develops the new imaging code based on a static 
point scatterer. Part 3 and 4 will further develop the code to cater for different waveforms 
and applications to real data. In this thesis, the physics involved and the approach for 
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addressing image artifacts common to targets moving in an unknown fashion will be 
discussed. The study of the new image model through correlation of the image data to the 
radar model will also be analyzed and developed via Matlab simulation. This image 
model will be compared with the corresponding point-spread function to exa mi ne the 
localization in phase-space (position-velocity). The behavior of the imaging model will 
also be examined under different aperture geometries. 

E. THESIS ORGANIZATION 

The thesis is organized in the following manner: 

Chapter II illustrates the development of radar imaging theory by looking at the 
basic building blocks necessary for imaging—radar operating principles, scattering 
models, correlation reception—and how these components translate to imagery. Finally, 
the thesis includes an overview of the inverse nature of the radar imaging problem to 
provide better appreciation of issues pertaining to practical measurement systems and 
how information is extracted from measurements to obtain the best estimates of the 
object. 

Chapter III applies the concepts put forth in Chapter II to the imaging scenario. 
The discussion here is based on the work of Cheney and Borden [2] in imaging moving 
targets from scattered waves. The effects of multiple moving targets on the scattering 
model are described. The approximations invoked, specifically slow-moving, narrow- 
band, and far-field, are examined. A scattering model to be used in the simulation 
program will be presented. 

Chapter IV will discuss the radar ambiguity function and the relationship between 
estimation errors and the transmitted waveform. Particular emphasis will be given to 
image quality analysis and the derivation of the point-spread function to quantify the 
performance of the imaging scheme developed in Chapter V. 

Chapter V exa mi nes further implications of the transmitted waveform through the 
pulse compression perspective. In particular, the Chirp signal is discussed. 
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Chapter VI demonstrates the behavior of the imaging scheme. A simple 
simulation program is developed and the results of the simulation, specifically the effect 
of different aperture geometries, are discussed. This chapter also summarizes the major 
findings of the thesis work and provides recommendations for future work. 
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II. IMAGING THEORY 


In Chapter II, we will look at some of the basics in imaging theory and examine 
how spatial images are formed. We begin by examining the physics behind extracting 
range and Doppler shift (range-rate) information using a radar system. The chapter 
concludes with a brief look at the inverse nature of the imaging problem and 
characteristics of measurement systems and some of the issues and challenges pertaining 
to radar imaging. 

A. RADAR SYSTEMS 

Radar systems can be classified by two main designs, namely stepped-frequency 
and pulsed. Most of the fielded radar systems are pulsed systems, as opposed to stepped- 
frequency systems. There are two key measurements that radar systems make: the 
strength of the backscattered field and the round-trip delay time of transmitted signal 
pulse reflected from distant targets. The range from the radar to the target is determined 
by measuring the time delay (round trip) for the radar pulse to travel to the reflecting 
target and back. Because the radar pulse travels at the speed of light we can deter mi ne the 
range R from R = ct I 2 where ris the round trip time delay. 

There are several aspects in understanding how radar signals provide information 
for imaging. Firstly, is the kind of waveforms that radar uses to transmit a signal. Radar 
signals are real-valued functions of time but it is convenient to express the waveforms as 
complex-valued functions. The radar system is able to measure both the amplitude and 
phase of these signals as functions of time and so a radar signal is of the form 
the function (p(t) is the “phase” and provides the phase information. The function A(t) is 
the “amplitude” and provides the magnitude of the waveform. 

All radar waves are electromagnetic waves governed by Maxwell’s equations and 
satisfy the wave equation and radiation condition. Essential target information such as 
range, speed, and bearing can be extracted from the waveforms. 
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The resolution in the down-range direction is determined by the pulse duration 
and the pulse bandwidth. A typical carrier frequency lies in bands which coincide with a 
low atmospheric attenuation window—l-40GHz. Larger bandwidth will imply finer 
resolution for imagery in the range direction. The concept of resolution will be addressed 
in Chapter III. 

B. SCATTERING OF ELECTROMAGNETIC WAVES 

In this section, we look at the behavior of electromagnetic waves and, in 
particular, the propagation and scattering of these waves. Knowing the properties of such 
waves provides important information about the transmitted and received signals and also 
provides target information. The following two sections will further develop these ideas. 

From the discussion in [18], a simple one-dimensional scattering model for a 
moving PEC^ plate in free space is a commonly used approximation for a typical radar 
echo wave. From Maxwell’s equations, it can be shown that an electromagnetic wave in 

the far-field consists of the electric field E , magnetic field H , and wave vector (direction 

A 

of propagation) k. These three vectors are mutually perpendicular. For a transmitted 

A A 

wave with its field component orthogonal to the x-direction with k = x, the boundary 
condition on the PEC results in the following relationship between the scattered and 
transmitted electric field: 

Escatt g{t + R{t)l C) =-Einc f(t-R(t)/ C) 


This scattered field consists of a left traveling wave (moving in the negative x -direction) 
where position of a phase-front is given by xit) = R(t). Consider a plate that moves 


linearly in the x-direction, so that R(t) = x + vt 


x=R(al boundary) 


. By making use of 


substitution M = t + R(t) I c can solve for t in terms of u at the boundary: 


^ A perfect electrical conductor (PEC) allows the charges to move freely and instantaneously in 
response to a field; consequently the fields inside a PEC are zero. The resulting boundary conditions are: 
the tangential components of the electric field must be zero and the tangential components of the magnetic 
field are related to currents flowing on the surface of the PEC. 
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u - Rl c 

t = - 

1 + v/c 


and obtain the expression for the scattered field: 


Escatt g(t + R{t) /c) = -Einc f[a(t + x/c-R/c)-R/c] (2.1) 


where a is the Doppler scale factor: 


1- V / c 

a =- 

1 + v/c 


( 2 . 2 ) 


For a transmitted field consisting of a signal waveform s(t) mixed with a carrier 
wave cosico^t), the scattered signal received at the antenna {x = 0) is: 

Pr^c^t) » sit -2RI c)cos [(1 - 2v / c)(t - Rl c)- Rl c]^ (2.3) 

This result is based on the assumption that the argument a(t - Rl c)- Rl c ~ t-lR! c . 
We also assumed that the scatterer is moving slowly so that / c is small. In addition, we 
can expand the denominator of Equation 2.2 in a geometric series and approximate a by 
1 - 2v / c. The multiplier of t in the argument of the cosine function is the carrier 
frequency (relative to the transmitted carrier frequency) shifted by an amount 

2v 

=-( 2 - 4 ) 

c 

This is the Doppler shift. As in the stationary scatterer case, the received signal is time- 
delayed, but there is also a Doppler-shifted version of the transmitted signal. 


C. CORRELATION RECEPTION-RADAR DATA MODEL 

One of the fundamental requirements for effective radar operation is detection of 
signals in noise. Radar detects the presence of an echo signal reflected from a target and 
extracts information about the target for imaging. The conservation of energy requires 
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that the wave intensity from a finite-sized antenna decreases in strength as the inverse 
square of the distance. As a result, the round trip signal power received by the radar 
decays by a factor of H relative to the signal that was transmitted. For example, at 
typical radar operating distances, the received power can be as low as 10'^^ Watts. This 
problem is further complicated by the presence of thermal noise at the receiver and may 
result in poor image formation. 

There are many ways to overcome this problem for example, increasing the power 
of the transmitter can help to increase the signal-to-noise ratio (SNR), but there are 
practical limits to relying solely on this approach. Coherent integration of multiple pulses 
is another method, but the main practical effect of an increased pulse repetition frequency 
decreases during the maximum unambiguous range, and the fact that targets motion 
during data collection and pulse integration alters the phase of the scattered field. Hence, 
the subsequent discussion examines the idea of correlation comparison performed 
coherently so that the phase of the transmitted signal is preserved in the reference signal. 

One way to allow coherent correlation is to transmit long coded pulses, together 
with appropriate signal processing techniques called matched filtering or correlation 
reception. This is known as pulse compression (see Chapter V). The correlation receiver 
or matched filter is an important example of a radar signal processor for detection of 
desired signals. The output of the matched filter is the cross-correlation function of the 
received signal and the expected received signal, corresponding to the transmitted signal. 
Hence, there will be a mathematical equivalence between the correlation receiver and the 
match filter receiver, and it is possible to implement the matched filter as a correlation 
process. 

Correlation reception introduces an intuitive approach for developing the standard 
radar data model used in imaging and the result is optimal detection in additive noise. For 
example, with a presence of Gaussian noise, correlation reception will select the best- 
guess estimate for the radar signal associated with the field reflected from a target, the 
scattered signal s^^Rt). From [19], we can see that such processing, when averaged over 

all time reduces to finding a function that minimizes: 
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(2.5) 


ikJ’i df+ df- 


2Re f 5 (t')s* (t')dt' 

j rec ^ ' scatt ^ ' 


The key to this approach is to determine from the random measurements 

of The first two terms are system-scenario specific. Searching for the best 

minimizing function hinges on the last term and can be made more efficient by restricting 
the search to a few parameters derived from some model. In this case, the natural model 
is based on the scattering interaction between the interrogating field and the target. If 
(t) denotes an incident pulse transmitted by the radar, then the linear radar (weak- 
scatterer) scattering model follows by superposition 

00 

= JJ P{o, T)s.Jt - T)e''^^‘-^^dTdo ( 2 . 6 ) 


where p{v,T) is the target reflectivity density for a radar system that measures two 
parameters, range and radial velocity. This reflectivity density is defined in such a way 
that p{o,T)dTdo is proportional to the field reflected from the target at range between 
CT / 2 and c(t + dr) I 2 with Doppler shift between o and o + do. 

From Equation 2.6, the scattered field is naturally modeled by the two parameters 
T and o. The search space will be parameterized by the time shift r and the frequency 
shift o . Considering the last term of Equation 2.5 and maximizing for all r, o in the 
scattered model, the correlation integral is: 

00 

?j(o,T) = \ T)e-^^^‘'-^^dt' (2.7) 


This correlation integral seeks out the component of the received signal that “matches” 
the time-delayed, frequency-shifted version of the transmitted signal. Unlike signals will 
be suppressed and, consequently, correlation receivers attempt to find the r and o that 
maximize the real part of Equation (2.7). The expected output of the correlation receivers 
is obtained by substituting Equation (2.6) into the correlation integral of Equation (2.7). 
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( 2 . 8 ) 


ri{v,T) = JJJ p{v\t')s.^P'- 


Notice that in Equation 2.8 the correlation noise term has been suppressed. The 
correlation noise term, which is the correlation integral between random noise and the 
time-delayed, frequency-shifted version of the transmitted signal measures how well the 
scattered signal matched the random noise will usually be small for coherent systems. 
Rearranging Equation 2.8, 


00 

; 7 (y,r) = Jj p(o\r')z(o-o\r- 

-00 


T')e ^ dr'do' 


(2.9) 


where 


1 1 

Z(o,t) = j s.Jt'--T)s*Jt'+-T)e"’‘'dt' (2.10) 

-00 

Equation 2.9 is the standard radar data model and it describes the output of the correlation 
receiver as the convolution of p and z up to a phase factor. The function ;i'(t>,r) 
defined by Equation 2.10 is the radar ambiguity function. Radar imaging is closely 
related to the radar ambiguity function that characterizes the accuracy of the target 
position and radial velocity from the radar data (see Chapter III). 

D. ONE-DIMENSIONAL (HIGH RANGE RESOLUTION) IMAGING 

One-dimensional imaging involves the generation of range profiles, which can be 
used to describe target substructure. Eor high range resolution (HRR) techniques, the 
transmitted pulse’s instantaneous range resolution is smaller than that of the target. As it 
sweeps across the target it sequentially excites the target’s scattering sub-elements which 
re-radiate energy back to the receiver. When these scattering sub-elements are non¬ 
interacting and point-like, the scattered pulse will be a sum of damped and blurred 
images of the incident pulse, which are shifted by time-delays that are proportional to the 
sub-element’s range [1]. Eigure 1 shows an example of one-dimensional images created 
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the HRR radar systems. It can be observed that the profile displays additional peaks 
outside of the target support, clear evidence of the problem associated with a single, non¬ 
interacting scatterer approximation. 


y (motors) 


0 in 20 :«) -10 



Figure 1. Example of a range profile from a B-727 jetliner. The top view (with 
orientation at time of measurement) is displayed beneath. From [1]. 


The signal received from a distribution of non-interacting, stationary targets can 
be modeled as 

00 

s^^p)=\ p{TXJt-T')dT'+n{t) (2.11) 

-00 

where n(t) represents an additive and random noise term. Application of a matched filter 
to obtain the radar data model used for imaging 
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00 00 00 

//(r) = j T)dt' = j z") j p(T')s.^^(t'- T')dT'dt' + noise term 

—00 —00 —00 

00 

= jj T)s.^^(t'-T’)dt'p(T’)dT' + noise term 

-00 

Making the substitution t" = t'-r' and neglecting the noise term (negligible with 
matched filtering) gives 

00 

r/(^)= J (2.12) 

-00 

where % is the autocorrelation 

00 00 

j(r)= j s\t"-T)s{t")dt" = \ s\t')s{t'+T)dt' (2.13) 

—00 —00 

Equation 2.12 represents a one-dimensional image and shows how the true distribution p 
is related to the radar data model. The image is the convolution of p with which in 
high-range resolution imaging is also commonly referred to as the point-spread function 
(PSF). The point spread function obtains its name from the fact that if p consists of a 

single point d{t), then the image is ?](t) = j xO -t')S(t')dt’ = xit ). Thus, x quantifies the 

degree to which a single point appears spread out in the HRR image. 

The down-range profile can be affected by a variety of factors such as the target 
aspect angle, position of the scatterers, or masking of scatterers by other parts of the 
target. Additionally, while the use of short pulses enhances the resolution, it concurrently 
leads to large bandwidth requirements. Wide bandwidths can increase system complexity 
and increase the likelihood of interference from other emitters in the electromagnetic 
spectrum. A short-pulse waveform also provides less accurate radial velocity 
measurement, a natural consequence from the properties of the radar ambiguity function 
(More details on radar resolution and the ambiguity function in Chapter III). An 
important implementation limitation in practical radar applications is the required high 
peak power to transmit short pulses over long ranges. High peak power transmission can 
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be problematic, especially at high frequencies due to the small waveguide dimensions 
and the small area of the anodes of the microwave sources. 

Radar target recognition using only range profiles has limited applications. This is 
because a range profile will not be able to distinguish cross-range target structures. All 
scatterers located at the same distance from the radar will reflect energy back to the radar 
with the same time delay. Hence, when the radar illuminates many distinct targets at any 
instant, meaningful interpretation of on-scene target (multiple) disposition cannot be 
derived based on a set (single) of range-only data. 

E. TWO-DIMENSIONAL IMAGING 

In order to interpret a target more effectively, additional information is needed on 
top of the range profile. This information can be in the form of a high-resolution cross¬ 
range profile, a Doppler profile, or simply the “triangulation” of range profiles. 

As a simple illustration for extension of the radar imaging concept to two 
dimensions, “triangulation” of different sets of range profiles will be explained (Figures 2 
and 3). This approach enables the determination of cross-range target structures while 
using only HRR radar systems and relies on collecting multiple sets of range profiles 
from different target orientations, processing them, and subsequently synthesizing a two- 
dimensional image. 

A simple example is used to illustrate this approach. Consider a set of three point 
targets with the radar located at the same distance from target 2 and 3 in Figure 2. When 
the geometry is oriented as shown, the return echo will only indicate two targets. Hence, 
ambiguity exists when targets lie along the bands of constant range from the radar. 
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Target 1 

o 


Bands of constant 
range from radar 


Target 3 


o 


o 


Target 2 


Single pulse ambiguity 


Figure 2. Ambiguous scenario from a single radar pulse. From[19]. 

Using multiple sets of data taken from different directions, triangulation will 
allow for gradual buildup of the relative positions of the three targets, as shown in Figure 
3. The range profiles are swept in the cross-range direction to form bands of constant 
range from the radar; these bands represent the possible locations of the target scatterers. 
The crossing points of these bands are used to determine the scattering center locations. 
The success of such an imaging scheme hinges fundamentally on the ability to correlate 
or correctly superimpose the various constant range bands. Otherwise, the intersection of 
the swept lines will be mistaken, resulting in image artifacts. 


Multiple pulse ‘'iuiagiiig” 


Target 1 


o 


Target 3 


c 


o 


Target 2 



Figure 3. Cross-range information obtained from range profiles. From [19]. 
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When a radar system is used for imaging a moving target, the received signal is a 
time-delayed and Doppler-shifted version of the transmitted signal. The measured echo 
signal from a single point-like moving target can be expressed as 

~ + ^^(0 (2.14) 

Because there are two unknown parameters to estimate, namely the delay r and 
the Doppler shift o , a single matched filter will no longer be enough. A set of matched 
filters, one for every possible Doppler shift, will have to be applied as 

00 00 

? 7 (r,t>) = j h^(t = /? J h^{t - -i- noise term 


in order to estimate r and o . Estimated values occur where I rj I takes on its maximum. 
For a distribution of non-interacting moving targets, the output from the filter bank is 

00 

rj(T,u) = j slit'- T)e-‘^‘-‘'-^^s^Jt’)dt' 

-00 

00 00 

= I II piT',o')s.^it'- do'dt' 

-00 -00 

00 

= III IJP- T')e-‘'’'^''-^'^dt'piT',v')dT'do' (2.15) 

-00 

00 

= III slit"+ t'- _^(t'Oe-'"'‘''Vt>(r',y')dr'dy' 

—00 

00 

= II Xi'^~ piz' ,v')dT' do' 

-00 


where 


ZiT,o) = I sl(t"+ T)e^"’^‘"^s(t")dt" = ^ 4e(^'+ (2.16) 

-00 -00 

is the radar ambiguity function. Similar to the one-dimensional case, the ambiguity 
function can be interpreted as an imaging point-spread function for range-Doppler 
imaging. Analysis of the ambiguity function determines the fidelity of the range-Doppler 
imaging process and consequently the accuracy of target range and velocity estimation. 
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Range-Doppler Spatial Image 


Range and range-rate measured from radar returns can be interpreted to form 
spatial images. The following example will illustrate this, provided the relative motion 
between the antenna and the target is known. Consider a simplified two-dimensional 
geometry between the radar and target as shown in Figure 4. The target is assumed to 
have translational and rotational motion only in the two-dimensional plane relative to a 
stationary radar platform. 

Consider a rotating target [18]. For simplification, assume that any overall 
translational motion is removed. 



lUuliir I’ulse 

l{(t) = I1(\+ (I siiif? 

▼ 

Figure 4. The geometry for locating a point target known to be rotating. From [18]. 

A point target rotates counter-clockwise in geometry in which the antenna is in 
the plane of rotation (Figure 4). The coordinate origin coincides with the center of 
rotation and the center of rotation is a large distance from the radar antenna. The large 
range approximation means that the incident wave is approximately a plane wave. 
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Consequently, all horizontal points are approximately at the same range. If the target is 
positioned at (x, y) = a(cos6',sin6'), then its range is 

R = +y = +asind (2.17) 

g g 

If its angular velocity is O-dOldt, then its velocity vector is a^(-sin^,cos6’), its 

g 

tangential speed is \v\=a6, and the down-range (y) component of the velocity is 

g g 

R=\v\co?,0 = aOcosO. However, since x = acos0, this down-range velocity can be 

g g 

written R = xO, where x is the cross-range component of the position. Note that for 
X > 0, the down-range velocity is negative, whereas for x < 0, it is positive. 

This means that because a radar system can measure range -i- y) and range- 

rate {R = x6), it can determine the (x,y) location of the point in space provided » is 
known. In particular, the radar measures the time delay r = 2R! c and the Doppler shift 

2R 2x6 

Vd= -^0=- 

c c 

which means that the coordinates of the point are given by 

(x,y) = (-^^, (2.18) 

^0 2 6 2 

This example shows how the range-rate radar measurements can be used to 
estimate Doppler and cross-range information. Further, the simple example highlights the 
difficulties in practical implementation of the radar scattering model. Estimates of the 

A 

scattered field require any variations in range R to be estimable. While it is reasonable to 

A 

assume in the one-dimensional imaging methods that R is constant over the data set, this 

A 

is not true for the two-dimensional imaging schemes. R will generally vary due to the 
target (translational) motion that occurs during measurements made while the target 
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rotates. The effect of failing to correctly account for this translational motion will be the 
introduction of phase errors into these data and a source of image corruption. 

F. RADAR IMAGING—AN INVERSE PROBLEM 

A discussion of imaging theory will not be complete without addressing the 
inverse nature of the imaging problem and its associated issues. As this thesis is focused 
primarily on developing and analyzing the point spread function rather than imaging, the 
intent of this section is to highlight the key characteristics of such inverse problems— 
issues pertaining to practical measurement systems and how information is typically 
extracted from measurements to obtain the best estimates of the target object. 

Radar scattering can be well-approximated by a linear, shift-invariant system. 
This is an important characteristic, as it states that the object function f{x) at each point 
X has a corresponding output image m{x). Further, the system will measure objects the 
same no matter when they start or, equivalently, where they are located. Radar imaging 
will not be meaningful if the delayed version of the target looks different for different 
delays. A general formulation of the imaging problem can be achieved by considering an 
imaging system represented by the functional operator or “kernel” k that describes how 
the measurement system works 

m = Kf (2.19) 

The direct problem in radar imaging refers to the mapping from the target to the 
quantities that can be measured by the radar. Consequently, the inverse problem is 
concerned with “reproducing” the original target from given data and knowledge of the 
direct problem; in the case of radar imaging, it will be radar measurements and 
knowledge of the scattering model. 

I. Well-Posed and Ill-Posed Problems 

One of the issues for inverse problems is their ill-posed nature. An ill-posed 
problem is one whose solution is not unique and/or does not exist for any data and/or 
does not depend continuously on the data. In the case of a band-limited system, the 
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solution of the inverse problem is not unique. This is because the imaging system does 
not transmit information about the target at frequencies outside the band of the measuring 
system. 

Consider discretized measurement systems in which the measurements are 
collected at specific points, and so Equation 2.19 can be written in matrix form. The basic 
goal of getting a good estimate / from the measurements m can be notionally 
represented as 

7 = K^m ( 2 . 20 ) 

It can be observed that k , which represents the measurement system, imposes certain 
properties on the measurement system. In general, the matrix k isA/^xM(M>A^)as 
the object function is infinitely dimensional and will not be confined to the limitations of 
the measurement system. Consequently, Equation 2.20 has more unknowns than there are 
linearly independent equations and so the system cannot have a unique solution. 

Since the object space (M-dimensional) is typically larger than the measurement 
space (N-dimensional), there exists what is known as the nullspace which consists of all 
vectors / - Kf or, physically, all the things that cannot be measured. Null space 
accounts for measurement artifacts and is determined by the kernel k . 

The accepted approach for solving inverse problems that are ill-posed is to search 
for approximate solutions satisfying additional constraints based on the physics of the 
problem [22]. This set of approximate solutions corresponding to the same data function 
is the set of objects with images close to the measured one and is expressed in what are 
known as the “normal equations” 

K^Kf=K^m ( 2 . 21 ) 


2. Data Reconstruction—Regularization 

Equation 2.21 specifies the set of conditions that must be satisfied by the least- 
squares solution / which can be rewritten as 
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/ = {K^KY^K'^m ( 2 . 22 ) 

Recall that the measured data suffers from finite dimension and noise contamination. 
These are problems since k) is poorly conditioned and the effect of noise is 

magnified (i.e. errors increase) as the dimensionality of the problem gets larger. Two 
techniques will be highlighted here, namely the singular value decomposition (SVD) and 
Tikhonov regularization, to mitigate the effects of noise when estimating / . 

From linear algebra methods, the matrix (k^kY^k^ can be expressed in terms of 
its SVD [19]: 

(k^kY'ic^ = UDV^ (2.23) 

where the matrices U, V contain the eigenvectors of , and D will have diagonal 
elements of the form 1/ (where A,, are the corresponding, non-zero eigenvalues). A 

consequence of the kernel k being bounded is that the eigenvalues necessarily form a 
sequence that gets arbitrarily close to zero as the number of measurements increase. The 
diagonal elements of D get large when higher dimensional data are used, and if the data 
are contaminated by noise, it follows from Equations 2.22 and 2.23 that noise will also be 
multiplied by a large value. Hence, one simple approach to mitigate the noise effect is to 
truncate the SVD representation when A. falls below some chosen threshold value. The 

corresponding diagonal element is set to zero. This method is known as the truncation 
filter. 

Tikhonov regularization modifies (/e^/e) to a matrix with better condition 
number by adding a fixed scalar a to the original normal equations. Equation 2.22 can 
then be rewritten in the form 

f = (K^K + aIY'i^^m (2.24) 

When a = 0, the original normal equations are returned and so it can be expected that 
small values of a will not change the system description too much. 
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After introducing a , the diagonal elements of D in Equation 2.23 become 



/I. + « 

l 


when 1. » a , D . « 1/ JA. (as before). When A. = a, however, D ~ JA. / 

^ ii \j i i ii i 

not become arbitrarily large. Consequently, the noise problem is moderated. 


(2.25) 

a and will 
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III. IMAGING DATA MODEL 


So far we have discussed the different types of imaging theories and techniques 
used for radar imaging. This chapter will make use of these ideas and concepts to develop 
an approach to imaging moving targets from scattered waves. We will also be looking at 
the three-dimensional case in order to develop the imaging model. 

A. INTEGRAL EQUATION APPROACH TO SCATTERING 

From the previous chapter, we have seen that for the one-dimensional scattering 
problem from a moving plate, the signal that returns to the radar is a time-delayed and 
frequency-shifted version of the transmitted signal. This is also the case for three- 
dimensional scattering from a moving, point-liked target. A point target is the traditional 
model used for developing imaging analysis and will also be used for the thesis. The 
representation of a target as points is appropriate because for high-frequency waveforms 
scattered from smooth conducting surfaces, the net scattering effects from the surface is 
equivalent to the contribution from the entire series of points known as specular points 
[19]. 

Applying traditional electromagnetic approaches to scattering will typically 
involve solving a differential wave equation and imposing the appropriate boundary 
conditions at the surface of the scatterer. However, for realistic radar targets, this process 
is difficult to implement because there is no coordinate system that can represent the 
boundary conditions for complex targets. 

Therefore it is convenient to reformulate the scattering problem. From the 
mathematical approach, it turns out that there is an equivalent integral equation for every 
differential equation, and this integral equation can inherently incorporate any boundary 
conditions in the solution. 
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This concept can be explained as follows: take a scattering object to be located in 
free space with an index of refraction n(x). Using a monochromatic wave, with harmonic 
time dependence to illuminate on a object. The reduced wave equation of the field 
becomes 

V^'¥(x) + k^n\x)W(x) = 0 (3.1) 

This is also known as the Helmholtz equation and k is the wave number in free space. 
Note also that n(x) on the scatterer. By splitting the field 'T into the incident field 

'T. which is the incoming wave, and the scattered field 'F , Equation 3.1 can be 
simplified to 

= P(x){^i„c(x) + ^,caAx)) 

where the scattering density p(x) = k^ (l-n^(v)) is a source factor completely 

determined by the index of refraction of the scatterer and is zero outside of the scatterer. 

The solution to Equation (3.2) can be obtained based on the application of Green 
functions, and we can write the solution as 

This is a Lippmann-Schwinger equation. It can be observed that the problem of solving 
Equation 1.1 has been reduced to finding the Green function G^{x',x) in free space (k 

denotes a specific frequency). The scattered field is also written as an integral over the 
scatterer D since p(x) = 0 otherwise, and hence automatically satisfies the required 
boundary conditions for a known target. 

It is also important to note that there are several aspects to the Lippmann- 
Schwinger equation: 
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(1) There are generally many functions that can serve as Green functions 
Gj^(x\x) in free space. While there are no bounding surfaces other than those defined by 
the scattering object, we are interested in G^(v',v) that obey the Sommerfeld radiation 

condition \\mGAx\x) oc - as this is more applicable to the typical radar operations. 

r^oo f* 

(2) The inverse scattering problem, i.e. that of determining p given and 
^, can also be approached using the Lippmann-Schwinger equation as a model for 

^ . 

scatt 


B. LINEARIZED DATA MODEL (TIME-VARYING SYSTEMS) 


Besides using the integral equation approach to solve the scattering problem, there 
is also a need to make the appropriate approximation in order for the radar model to be 
tractable. By enforcing the Bom or “weak scatterer” approximation, where the object 
distribution (cell) is not self-interacting, multiple scattering events are treated as 
negligible in comparison with the primary scattering events. The result of this is a 
linearized (approximate) data model as follows: 

(3.4) 


Equation 3.4 turns out to be a solution for Helmholtz’s equation. A typical radar 
operation is a time-varying problem, and so the time-dependent wave equation and the 
time-domain Green function will have to be incorporated with the model 


Gi^(x\x)^ g(x\x;t\t) 


S(t -t'- I x'- X I /c) 
An I x'- X I 


(3.5) 


With a signal s^(t) transmitted from position y at starting time -T^ , the incident field at x 
is 


5 (t -r T - \ x- y\lc) 

y/. {x,t)--— -- 

4;rlx-yl 


(3.6) 
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Substituting the incident field (Equation 3.6) and the time-domain Green function 
(Equation 3.5) to Equation 3.4, we get the scattered field at time t and 

position z 




.(j,z,0 = jj 


d{t-t'- \ z-x' \ Ic) Sy{t'+T^- \ x'- y \ Ic) 


4;r I z - .r' I 


An\x'-y\ 


-p{x')d^x'dt' (3.7) 


m 

where Sy denotes the second time derivative of incident signal 5^.(0 and enters as a 

consequence of the second time derivative in the wave equation. Hence, as a result of the 
Born approximation and the usage of time-domain Green function, we obtain a linearized 
data model (Equation 3.7) 

C. REFLECTIVITY FUNCTION FOR MOVING TARGETS 


Recall that the reflectivity functions p{x ), or scattering density function first 
introduced in Equation 3.2 is determined by the target. It is a scale factor of the received 
signal strength and an important characteristic of radar imaging. Eor a moving target, a 
time-varying reflectivity function p(x,t) will be used. When the imaging scene contains 
multiple moving targets or scatterers, the scattering model can be modified as follows: 

Let p^(x-vt)d^xd\he the scatterers in the volume element d^xd\of phase 
space centered at position v and velocity v. Choose T , the time when the signal is first 

send out, to be such that the scattering density p^ = p^(x). This means that all the targets 
in the scene moves with velocity v at time t=0 (i.e, the signal is transmitted from the 
transmitter antenna at time t = -T ^). As a result, the spatial scatterer density centered at 

time t and position x is 

p{x) = ^ p^ix - vt)d\ (3.8) 

and the scattered field becomes 
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Vs, 




5{t-t'-\z-x'\lc) Sy{t'+T^-\x'-y\lc) 


A7t\z-x'\ 


An\x'-y\ 


pXx'-vt')d\d\'dt' (3.9) 


In order to simplify the problem, making a change in variables x' ^ .r = a:'- vt' can 
change the frame of reference in which the scatterer is fixed. As a result the scattered 
field is modified to 


Vs, 




5{t-t'-\x + vt'-z\lc) Sy{t'+T^-\x + vt'-y\lc) 


4;r I V + vtz I 


4;r I X + vty I 


pXx)d\d\dt' (3.10) 


The data model, or scattered field, is now modified to incorporate moving targets. The 
physical interpretation to Equation 3.10 is that the wave that emanates from radar 
transmitter position y at time -T encounters a target at timet'. This target, during the 

interval[0,t'], has moved from v tox + vt'. The wave scatters with strength p^(x) and 
then propagates from position x + vt' to z, arriving at time t. 

To further simplify the equation, let 


R^^(t) = x + vt-z, R=\R\, R = R/R 


Substituting this to the scattered field Equation 3.10 gives 


Vs. 




S(t-t'-R^(n/c) Sy(t'+T-R^(t')lc) 


ATzRid) 


ATzRid) 


p^(x)d^vd^ xdt' 


(3.11) 


The variable t’ in the argument of the delta function appears in several terms making it 
difficult to carry out the integral over t’. To address this, let t' = T^ ^(t) denote the implicit 

solution of 


t-f-R^^(f)lc = 0, i.e., t-t^ ^(t)-R^ ^(t^ ^(t))l c = 0 

Then, performing the integration over t' yields 


29 



V^scati(y’^’0 = -^-=-=- -p^(x)d vd X 


Substituting t^^{t) = t - ^(t^^{t))! c for the first occurrence of in this last 

equation allows for a more symmetric version 


¥scaniy^Zd) = \\ 




p^{x)d^vd^x 


(3.12) 


D. FURTHER SIMPLIFYING APPROXIMATIONS 

Equation 3.12 applies to very general scenarios. It applies to rapidly moving 
targets without the need for start-stop approximation, and it is appropriate for any 
transmitted waveforms. 

However, this model still has “retarded-time” problems. These problems arise as a 
consequence of the definition of retarded time which is the delay between the 

scattered fields reflected from the target to the radar receiver. The retarded time is 
expressed as t^^{t) = t-R^^(i^^{t))! c where t is the time when the scattered field 

reached the radar receiver. Hence, a more accurate result can be obtained when the 
scatterers are known to be moving in such a way that R^^it) can be expanded in a Taylor 

series around t=0 and approximated by retaining only the terms linear in t — i.e., when 
the scatterer is “slow moving”. 

I. Slow-Mover Approximation 

Assume that (I v 11) and {\v? f x co I c) are much less than \x- z\ and 
IX - y I, where ct>^^ denotes the maximum angular frequency of the transmitted 
signal 5^. In this case, 

R^,(t) =\z-ix-vt) l=7?^^(0)-i-Rx,z(0)g^t-f... 
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where R^^(0) = x-z, i?^^(0)=l/?^^(0)l, and Rx,z(0) = R^^(0)/R^^(0) . Substituting 
this result into the definition of retarded time yields 


t (0 « t - R,, (0) + Rx,z (0)q^K^ (0 


/C : 


t-Rx,z(0)/c 

1 + Rx,z(0)g^ /c 


Inserting this approximation into the result for y/ (y,z,t) yields 


¥scaAy^Z,t)=\\ - .. Vd X 


(47rrR,J0)R^J0) 


(3.13) 


where 


a = 

x,v 


l-Rx,y(0)g^ / C 
l + Rx,z(0)g^/c 


(3.14) 


is the Doppler scale factor. Similar to the traditional Doppler shift, it is closely related to 
the velocity component in target-transmitter direction. 


2. Slow-Mover and Narrow-Band Approximation 

Typical systems are narrow-band for which the transmitted signal is of the form 


where Sy(t) is slowly varying (as a function of t) in comparison with e , where o)^, is 
the carrier frequency for the transmitter at position y. In this case, the time derivatives of 

gg 

5^ are dominated by the e factor, and S y(t) » -co^ S y(t)e^““- '. Moreover, because 
Sy(t) is slowly varying, further approximations can be made 


Sy{a^J)e^ 


Sy(t)e 


Equation 3.13, under the slow-mover and narrow-band then becomes 


(3.15) 
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(3.16) 


CC UJ e e 

WscaAy^ Z,t) = -jj 

xSy(t + T^- (0) + (0)) / c) p, (x)dVx 

where 



(3.17) 


3. Slow-Mover, Narrow-Band and Far-Field Approximation 

When the transmitter-to-target and target-to-receiver distances are large in 
comparison with the scene dimensions, then l.r + vtl and \x + vt? xa> Ic can be 
assumed to be much less than either I z I or I j I. And the expansion 

A 

Rx,z =\z-(x- Vt) 1=1 z I - Zg;x + vt) + ••. 

can be applied and similarly for R^^{t). Substituting the expansions into previous 
approximations Equation (3.16) then yields 


Wscattiy^Zd) 


-CO„ 


(4;r) I z II y I 


11 


i(p, „ J 

e e 


(3.18) 


xS, 


t + r - (I z I - zgr+1 y I - ygr) / c p^(x)d^vd^x 


where now 


(D =(o 

T X.V V 


lyl-ygr-cr +«,(lzl-zgr) 


! c 


(3.19) 


Note that in Equation 3.16, the expansions for R^ , R^ ^(t) are used differently. In the 
amplitude factor, the approximation R^ =\ z - (x + vt) \»\ z \ holds since 

A A 

I z I? I z^x + vt) I; in the phase factor, however, the product (zgr x / c) can still be a 
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large fraction of 2;r, and must be retained. In addition, the signal envelope S can still 

A 

vary significantly over I z I - zgv despite being slowly varying, hence the first order term 
is retained. 

The Doppler scale factor also reduces to 

A 

1— yw I c A A 

«1 - ( z)g’ / c (3.20) 

1 + Zg^/c 

where the binomial expansion is used since I v I /c = 1. The quantity x is the 

A A 

Doppler shift (where = -(y+ z)g^ /c ) and is observed to be dependent on the “bistatic” 

A A 

vector (y+ z). 

Inserting =1 y I /c and setting = (0^1 c 


(p^^=k\z\-k(y+z)ix + {zi^z-x))v I c 


(3.21) 


E. IMAGING VIA A FILTERED ADJOINT 

In developing the correlation receiver, it was observed that the “fit” between two 
complex-valued functions is given by the cross correlation. For correlation reception, the 
correlation integral seeks out the component of received signal that “matches” the time- 
delayed, frequency-shifted version of the transmitted signal. Large values of the 
correlation integral indicate a strong resemblance, while small values indicate weak 
resemblance. The values of the time-delay and Doppler shift parameters, which 
maximized the correlation integral, were also the parameters most likely to represent the 
target. 

Using the same notion, imaging can be achieved by determining the position and 
velocity parameters that represent the best “fit” to the measured data. These will also be 
the parameters that localize the position and velocity of the unknown scatterers in phase 

space (position-velocity) space. Cross correlation in the time domain, is multiplication in 
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the Fourier transform domain (Fourier correlation theorem), hence the association with 
filtering since filters are multiplication operators in the Fourier transform domain. 

The model developed for y^,^^„(y,z,t) depends on the position of the scatterer 
X = (x^,x^,x^) and its velocity v = (v^,v^,v^). Using this form for y/^^^^^(y,z,t) to build a 
six-parameter scattering model based on the arbitrary position p = (p^,p^,p^) and 
arbitrary velocity u = (u^,u^,u^), the objective will be to determine the values of p and u 
that maximize the cross correlation between the arbitrary model and the measured data. 


In terms of the Doppler shift, the scattered field from known scatterers p^(x) with 
positions x and velocities v is 


V^sca„(y’Z,t)^ 


(4;rflzllyl 


JJ exp (y+ z)ix -v(t-zi^z-x)! c)] 


x5, 


t-(lzl-(y+z)gr)/c p^(x)d^vd^X 


The effect of data shifted in time is an image that is translated in range. Because 
this has no bearing on the imaging process, the scattered field result can be simplified by 
substituting t' = t-\z\lc which yields 


¥ scanty 


2 -io>J 

-G)ye 

(47r)^\z\\y\' 


jjexp -ik^Xy+ z)ix-v(t T zgr/c)] 


x5, 


y 


(3.22) 


t'+(y+ z)^ / c jp^(x)d vd X 


Ignoring the intensity prefactor, a parametric model representing the field from an 
unknown point source is 


y/(y,z,t')=-e ^ txpl-ik (y+z)ip-u(t'+z^!c)] 


(3.23) 


x5, 


t'+(y+z)&lc 


And the image is created as follows: For each p and u, determine the function 
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l{p,u) = jjj i//^^^^^{y,z,t')y/\y,z,t')dt'd’"yd"z 


(3.24) 


where m, n depends on the configuration of the transmitter(s)/receiver(s), and the 
integrals are over all values of t' ,y and z for which measured data is available. 

I(p,u) = jjj y/ ^^^^^{y,z,t')y/\y,z,t')Q{(o,t\p,u,y,zW d"‘yd"z 

. ,, iK(y+z)i p-u(t'+zgp/c) 

(J 0 ,t ,p,u,y,z)e e 

(3.25) 

t'+(y+z)s}/c]y/^^^„(y,z,t')dt'd'"yd"z 
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IV. IMAGING POINT SPREAD FUNCTION 


This chapter begins by looking at the radar ambiguity function its relationship to 
the transmitted waveform, and the accuracy with which target range and velocity can be 
estimated. The imaging point-spread function (PSF) is also derived based on the new 
imaging scheme. 

A. RADAR AMBIGUITY FUNCTION 


For a stationary target, the output of the correlation receiver is the cross 
correlation between (1) the received signal plus noise and (2) the transmitted signal. In 
many radar applications, however, the target is moving and so its echo signal has a 
Doppler frequency shift. The output is therefore a cross correlation between the received 
signal and the Doppler-shifted transmitted signal. 

The nature of the correlation receiver output as a function of both time and 
Doppler frequency is important for understanding the properties of a radar waveform. In 
particular, the choice of waveform affects measurement accuracy, target resolution, and 
ambiguities in range and radial velocity. These aspects of the correlation receiver output 
will be examined in turn. 


When the received echo signal is large compared to noise, the output of the 
correlation receiver may be written as (see Equations 2.15 and 2.16) 


;7(y,r) = JJ p(u',T')z(o-u',T-r')e 


i—(u+u')(T-T‘) 

^ dr'dv' 


where - o',t - r') is the radar ambiguity function. The problem of radar imaging is 
concerned with estimating the location and strength of the point scatterers that are 
assumed to make up the target. Consequently, understanding x is essential as the end 
state involves “inverting” the radar data model. By definition, we have 


X{v - v',x- t') 


1 






r. 1 


li'. I 

I me I 


t'+ —(r 
2 





dt' 


(4.1) 


-00 
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Basic Properties 


For simplicity, the origin is chosen to be the true target time-delay and frequency- 
shift. This assumption will be carried throughout the rest of the discussion. Some key 
properties that follow from Equation 4.1 [Refs. 14, 18 and 19] are listed 


Signal energy : \x{o,t)\< ;r(0,0) 


(4.2) 

Ambiguity volume : JJ ;jf(t;,r)| drdo^l 

-00 

{normalized signal) 

(4.3) 

Symmetry: \x{-v-r)\ = \x{v,T)\ 


(4.4) 

1 ” 

Frequency - domain: xi^^ = ~ 

In ^ 

—1>)5'* {co + —u)e '‘“''dco 

2 2 

(4.5) 


Equation 4.2 states that the maximum value of the ambiguity function occurs at the origin 
(true target time-delay and frequency-shift). Equation 4.3 asserts that the total volume 
under the ambiguity surface is a constant. This gives rise to a radar uncertainty principle: 
choosing a signal sit) so that the ambiguity surface will be narrow in one dimension will 
cause it to be correspondingly wide in the other dimension. In the case of range-Doppler 
imaging, a waveform with good range resolution has poor Doppler resolution and vice 
versa. 


2. Range Resolution 


Equation 4.5 presents an intuitive look at range resolution. Consider a fixed target 
(y = 0) and rewrite Equation 4.5 as 


I j(0,r) 1=;^ 
Ln 


fl5. 

J me ^ ^ 


(4.6) 


Recognize the form of Equation 4.6 as the inverse Eourier transform of the transmitted 
power spectral density. Then, in order for I ;j'(0, r) I to have high resolution in r 
measurements (ideally a delta function), the power spectral density has to be identically 
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one; such a signal will be broadly supported in the Fourier domain. Hence, it is clear that 
better range resolution is obtained when the signal bandwidth is large. 

A more rigorous study of estimation errors in the presence of Gaussian noise is 
developed in [1] by examining the behavior of the ambiguity function in the 
neighborhood of its main peak. This study shows that estimation error of o and r is 
described by an ellipse, and that time-domain resolution is inversely proportional to 
frequency-domain bandwidth. Hence, bandwidth is inherently important to problems in 
radar target imaging. A further interesting result makes uses of the definition of 
bandwidth in terms of the second moments (probabilistic variance) of the power spectrum 
for a narrow-band signal s(t) = . We can write 



da(t') 

dt' 


dt' -f 


1 

2^1 dt' ) 


\a{t') P dt'-al 


(4.7) 


where fd represents bandwidth. This useful relationship illustrates that a nonlinear 0(t) 
(phase modulation) will increase the bandwidth in comparison with signals that are 
modulated in amplitude only. This is the basis for pulse compression (see Chapter V), a 
technique of practical importance because it allows for the creation of fine range- 
resolution waveforms that are also of long duration with sufficient energy on target. 


3. Doppler Resolution 


To study Doppler resolution, consider a target whose range is known. Thus, 
setting r = 0 in Equation 4.1 yields 


I j(y,0) 1= 


□u 

fl5. (T)l 

J me ^ ^ 


2 iut 


dt' 


(4.8) 


It is evident that better Doppler resolution is obtained from a long duration signal. This is 
in stark contrast to range resolution requirements. In addition, the Doppler resolution is 
determined only by the amplitude modulation of the signal and not by the phase. 
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B. IMAGE ANALYSIS 


From Chapter III, an explicit form for the image is 


I(p,u) = -jjj Q{co,t\p,u,y, e 


, ik {y+z)i p-uit'+zgpic) 


xSy t’+(y+z)gp/c y/^^^,,{y,z,t')dt'd'"yd''z 


To examine the performance of the imaging scheme, the ideal (expected) data are 
inserted back into the image equation. Substituting for y/ (y,z,t') yields 


. rrrrrOltQiaj'.p.u.y.z) ^ ^ 

r A A A A 

X exp I ik^Xy+ p-x-(u-v)t '-u(zff>) / c + v(zgr) / c] 


(4.10) 


xSy t'+(y+ z)gp / c Si) t'+(y+ z)g!c / c d^vd^xdt'd'"yd"z 


The distances I z I and I y I are scaling amplitude factors and not essential for 
imaging purposes. Thus, choose 


Q{co,t',p,u,y,z) = 


(4;z-)" I z II y I 


J{p,u,y,z) 


(4.11) 


where J{p,u,y,z) depends on the geometry and is chosen to compensate for a Jacobian 
that results from the integral of the variable t’ 


• Mf. f>f> ~ ( A A \ ” f A A ^ 

/(p,m) = JJJJJsJ C+(y+z)gp/c S J t'+(y+ z)g>c/c \e 


-ik (y+z)g.u-v)t' 


xexp ik (y+z)^p-x-M(zgp)/c + v(zgr)/c] 


(4.12) 


X J(p,u, y,z)p^(x)dt'd^d^xd^yd^z 


Equation 4.12 can be written as 


I(p,u) = jj K(p,u,y,z)p^(x)d\d^ 


(4.13) 
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where 


K{p,u,y,z )=\\exp -iky{y+ z)^M(zgp)-v(zgr)]/( 


C " f AA \ ( A A 

JSj, t'+(y+z)s}/c Sy t'+(y+z)sc/c 


^-ikpy+z)Au-v)t' , 


(4.13) 


xj(p,u,y,z)d"‘yd''z 


is the point-spread function describing the behavior of the imaging system. Equation 4.12 
caters for multiple transmitters and receivers where y and z represent the transmitter and 
receiver location respectively. The corresponding superscript m and n accounts for the 
total number of transmitter and receivers. Apply the change of variables 

A A A A 

t = t'+ — (y+ z)^p + x)/c and set T = (y+z)^p-x)/c to obtain 


K(p,u,y,z) = \\ exp j -iky(y+ z)^ M(zgp) - v(zgr)] / c 


00 * 

r J, 1 

r 1 ^ 


( 1 ^ 

[ ^y\ 


Sy\ 

t — r 

J 1 

_-00 

1 2 ; 

1 2 j 


I A A ^ A A j 

X exp <1 ik^ iy+z)i{p-x) + -{u-v)iy+z)ip + x)]l c\ (4.14) 

The integral in square brackets is the radar ambiguity function 

e-‘^‘dt 

X is determined by the transmitted signal and I x I attains its maximum value at p = x, 
u = v for general configurations of transmitters/receivers. So the imaging PSF can be 
written as 


f ' Y 1 ^ 


r 

1 ^ 

[ Sy t + — T 


t- 

— r 

^ \ 2 J 

—00 X / 



2 J 


^(p, M, y, z) = jj exp j(y-i-z)^M(ygp) - v(zgr)] / c 


I A A 'Y A A j 

X exp I ik^Xy+ z)iip -x) + -{u- v)(y-f z)g;p + v)] / c 


(4.15) 


X J 


kSy+ z)^u - v),{y+ z)i^p -x)!c J(p,u,y,z)d'"yd’'z 
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V. PULSE COMPRESSION 


Detection of signals in noise and resolving capability are both fundamental to 
radar imaging. As discussed in Chapter II and IV, the use of correlation reception and a 
transmitted waveform with nonlinear phase modulation provides means to enhance 
performance in detection and resolution. Pulse compression is a technique that fuses 
these concepts and involves transmitting long coded pulses, together with signal 
processing methods such as correlation reception. This is of practical importance 
because it allows for the creation of fine range-resolution waveforms that are also of long 
duration. This increases radiation energy directed on target without relying solely on 
boosting radar transmitted power. This chapter gives a physical interpretation of pulse 
compression and looks at a specific phase-coded signal known as the chirp. 

A. REVISITING CORRELATION RECEPTION 

Recall that in Chapter II we showed that the output of the correlation receiver is a 
correlation integral between the received signal and a time-delayed, Doppler-shifted 
version of the transmitted signal. It is illustrative to explain the phenomenon of pulse 
compression with a graphical example. Using the simple one-dimensional imaging case 
in Equation 2.12, the correlation receiver output can be written as 

00 

77(t)= \ s\tXJt'+t)dt' (5.1) 

-00 

which is a correlation between s and s . If 5 = 5 (reflectivity function is 1), Equation 

5.1 is called an autocorrelation. In Eigure 5. the signal at the top represents the 
transmitted waveform. The lower graphs represent the received signal for different shifts, 
according to Equation 5.1. When a signal shifted by t is multiplied by the waveform at 
the top and the product is integrated over f, the resulting number is plotted at position t 
on the graph to the right. 
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Figure 5. Phenomenon of pulse compression where the energy is concentrated at a 

single delay time. From [18]. 

Figure 5 shows that although the output signal has a longer time duration than the 
original signal s, its energy is more concentrated at a single delay time. It is a 
consequence of the Cauchy-Schwarz inequality that the highest peak occurs when the 
signals are not shifted relative to each other. The smaller peaks, called sidelobes, are 
undesirable and an open problem is to design pulses whose autocorrelations have the 
lowest and fewest sidelobes. 

B. PHASE CODING—CHIRPS 

A radar system makes the most efficient use of power when it transmits a 
constant-amplitude waveform. Consequently, most radar waveforms are designed with 
variations only in the phase of the signal. The most important and most commonly-used 
radar waveform is the chirp. A chirp is a constant-amplitude signal whose instantaneous 
frequency varies linearly with time; chirps are therefore also called Linearly Frequency 
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Modulated (LFM) waveforms. Linear variation of the instantaneous frequency implies 
that dcj)! dt = + yt. The coefficient y is called the angular chirp rate. A chirp is thus 

of the form 


1 7 

, . ‘(“min' + Tr' ) . , 

s{t) = e 


(5.2) 


where j-](0 is 1 in the interval [0, 7] and 0 otherwise. 

A chirp with positive chirp slope is called an upchirp; one with a negative chirp 
slope is a downchirp. The instantaneous (angular) frequency of the chirp varies from 
CO . to CO = CO . + yt, with a center frequency of co,=co. +yt I 2. Thus the 

instantaneous frequency of the chirp varies over the interval - 7/ / 2,(»g + 7/ / 2] . 
The power spectrum of a chirp is roughly constant over the frequency band 
[®q-7//2,(»jj + 7//2], thus the angular frequency bandwidth Ty can be estimated 
by looking at the range of instantaneous frequencies (see Figure 6). 



H 11 
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Figure 6. Spectrum for a chirp pulse of length r . From [1]. 


Chirps are very commonly used in radar imaging. Pulse compression for chirps 
has the following intuitive interpretation: When the pulse is transmitted, different parts of 
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the pulse are coded with different instantaneous frequencies. In the pulse compression 
process, the different frequencies are delayed by different amounts of time, so that all the 
energy emerges at the same time [18]. 

C. AMBIGUITY FUNCTION FOR COMMON WAVEFORMS 


From the preceding discussion, the “ideal” ambiguity function will typically 
consist of a single peak of infinitesimal thickness at the origin and be zero everywhere 
else. It will closely approximate an impulse function, and have no ambiguities in range or 
Doppler frequency shift. The infinitesimal thickness at the origin will permit the time 
delay and/or frequency shift to be determined simultaneously to as high a degree of 
accuracy as desired. Such a highly desirable ambiguity function, however, is not 
theoretically allowed. The properties of the ambiguity function implicitly mandate that 
trade-offs exist between the measured dimensions. 

In practice, other than the trade-offs in dimension resolution, waveforms also 
generally have significant responses outside the narrow region in the near vicinity of the 
origin, i.e. there are many possible ambiguous (blind) range and Doppler. Ambiguities 
are a consequence of discontinuous waveforms, such as in a pulse train. The resolution 
for single rectangular pulse, single chirp and coherent pulse trains will be discussed next. 


I. Resolution for a Single Rectangular Pulse 


Time domain impulses have infinite bandwidth; hence they are idealized signals. 
More realistic waveforms can be modeled, for example, as 


Jl ifO<t<T, 

[O otherwise. 


(5.3) 


for which s. (t-\llT)s* (t + Hlr) is non-zero in the interval t e(lrl/2, r-lrl/2). 

me ^ ^ me ^ ^ V 7 / 

The ambiguity function is then 
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(5.4) 


f 1 . 1 ■ fT-Mn . ,, 

X{o,t)= [ 5 . it'--T)s if+-T)e"’‘dt'= [ e‘^‘dt 


= —e 
iu 


iu(T/2) 


^iv(T/2-M/2) _ ^-w(TI2-MI2) 


= {T-\t \)e‘^'^%mc 


-y(r-lrl) 

2 


For v = Q , then the delay resolution (by the peak-to-first-null measure) as 
observed from the ambiguity surface lj(0,r)P, is At = T (zero-crossing at t = T), 
which corresponds to a range resolution of AR = cT 12. For r = 0, the ambiguity surface 
I j(t),0) P yields Ao = llT because the first zero of the sine function occurs when its 
argument is n. The Doppler shift is related to down-range relative velocity by 
- -2vl)q / c (see Equation 2.4), and the corresponding velocity resolution is 
Av = IT. 


2. Resolution for a Single Chirp 


A chirped signal provides another example of a ridge-type ambiguity function 
produced by linearly frequency modulating a rectangular pulse over a bandwidth (5. The 
pulse width T is large compared to \l fd. The frequency modulation increases the 
spectral bandwidth of the pulse so that pT ? 1. Because the pulse width T and the 
bandwidth jd can be chosen independent of one another, the time-delay and frequency 
accuracies are independent of the other. 


The chirp signal is one whose frequency changes linearly with time 


5 . (0 = <; 

me ^ ' 


( .1 


ifO<t<T, 

otherwise. 


(5.5) 


The constant y is known as the chirp rate. The ambiguity function is then 
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00 1 

j(t),r)= [5. (?'— t)s* (t'+— T)e"'‘dt'= { 

' J me ^ y ^ me ^ ry ' J|j.|/2 


dV 


1 


i{v + yz) 
= {T-\T\)e 


iiu+r^)TI2 


^iiv+YT){TI2 -\t\12) _ ^-i{v+yT)iTI2 -\t\I2) 


(5.6) 




-{o + yz){T-\T\) 


Observe that a chirp leads to an ambiguity function that is a 5 /nc-function in both the time 
and frequency domain. It is instructive to examine this property again by looking at the 
signal resolution. 

For r = 0, a similar result of Ay = 1 / T is obtained as in the case of the 
rectangular pulse. This is because the amplitude modulation is identical in both cases. 
Consequently, Doppler resolution is the same. However, with phase modulation, it is 
expected that the delay resolution, and therefore resolution will be enhanced. 

For y = 0, when r>0 the first null occurs when n = (J -z)yzl 2. Reference 
[19] shows that the expression n = (T -z)yT 12 has a quadratic equation form with a 
solution l^z = 2n ! yT, which corresponds to a range resolution of lsR = cl /3 where 
P = 27zyT is the bandwidth in Hertz. For r = 0, Ay =1 ;t'(y,0) P= 1 / T because the first 
zero of the sine function occurs when its argument is tz . The Doppler shift is related to 
down-range relative velocity by y^ = “^vy^ / c (see Equation 2.4), and the 
corresponding velocity resolution Av = 2T . 

The pulse compression ratio is defined as the product of pulse spectral bandwidth 
P and the uncompressed pulsewidth T. Comparing the range resolution between the 
single rectangular pulse and single chirp. 


AR 


rectangular 

AR, 

chirp 


cT 12 _pT 
dp ~ 2 


(5.7) 


In other words, phase modulation improves the range resolution by a factor of one half 
the time-bandwidth product. 
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3. 


Resolution for Coherent Pulse Trains 


Following the analysis given in [18], a train of N identical pulses is considered: 


s(t) = 


1 


N-l 

^u(t-nT^) 


n=0 


(5.8) 


where is the pulse repetition interval, which might be on the order of 1ms. The 
reciprocal 1 / 7)^ is called the pulse repetition frequency (PRF). 

The ambiguity function for the pulse train is 


1 N-lN-l 

Z(o,t) = — X X j“*(^ + mT^)u(T'- nT^)e“'"'dT' 

^ n=0 m=0 


1 

N 


N-lN-l 

m)T^ , 


«=0 m =0 


(5.9) 


where Xu denotes the ambiguity function for u and where we have used the substitution 
t” = t'- nT^ . A lengthy calculation involving rearrangement of the order of summation 
and the summing of a geometric series results in 


Xiu,T) = \- 
N 




iuT (N-l-\p\) 


P=-(N-1) 


sin[;rL)r^(A- I p I) 
sin( ttuT) 


(5.10) 


If the pulses u are sufficiently well-separated so their ambiguity functions Xu do not 
overlap, then we take the absolute value of Equation 5.13 to obtain 


N-l 


N ' 

p=-(N-l) 


sm 


[TTUT^iN- I p I)] 


sin(;rt>r ) 


(5.11) 


Such an ambiguity function has a “bed of nails” appearance (see Figure 7), with peaks at 
T = pT^, p =-(A - 1),-(A - 2),...,0,1,...,(A -1) and u such that 

ttoT^N = (m + l I 2 ) 71 :, m = 0,1,2,... 
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Figure 7. Ambiguity Function for a train of pulses. From [2] 


Following the previous analysis, consider the u = 0 cut of the ambiguity surface 
to estimate the range resolution yields 

\z(0,t)\=^ ii \zM-pT^,0)\{N-\p\) (5.12) 


where we have used the fact that lim[sinZ?v] I sinx = b. We see that the range resolution 

x^O 

of the main peak is the same as the range resolution of the individual pulse u, but that we 
now have range ambiguities due to the extraneous peaks at t = pT^ . 


To estimate the Doppler resolution, consider the r = 0 cut of Equation 5.11 and 
for well-separated pulses, only the p=0 term contributes 


l;r(y,0)l=^ X 


sin[t;r^A] 

sin[yr^] 


(5.13) 


Consider the numerator argument because it varies much faster. It shows that the Doppler 
resolution is 1 / , and hence the Doppler resolution of the pulse train is higher than 

that of a single pulse. 
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The many ambiguities produced by a pulse train may appear to lead to a poor 
radar waveform. The time-delay measurement accuracy, determined by pulsewidth T and 
the frequency accuracy (determined by the pulse repetition interval ) can be selected 

independently. If pulse repetition interval is such that no radar echoes are expected 
with a time delay greater than , and no Doppler-frequency shifts are expected greater 
than 1 / , then the effective ambiguity function reduces to just a single spike at the 

origin whose dimensions are determined by T and . In practice, ambiguities can also be 
resolved with different PRF, and many radars employ this type of waveform. 
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VI. RESULTS AND CONCLUSION 


This chapter shows how the new imaging algorithm (Equation 3.25) is 
implemented and analyzed. Recall from Chapter III that the image is formed when the 
correlation receiver observed a “fit” between the scattering model (Equation 3.22) and 
the parametric radar model (Equation 3.23). Eor correlation reception, the correlation 
integral seeks out the component of the received signal that “matches” the time-delayed, 
frequency-shifted version of the transmitted signal. Earge values of the correlation 
integral indicate a strong resemblance while small values indicate a weak resemblance. 
The values of the time-delay and Doppler shift parameters that maximized the correlation 
integral were also the parameters most likely to represent the target. Using the same 
notion, three sets of Matlab programs were used for the implementation: 

1. The Scattering model (Equation 3.22) 

Scattering Model V^,catt^y,z,t) is basically the measured data from the scatterer. It 
provides all the information of the scatterer position x = (x^,x^,x^) and velocity 
V = . A code is built (refer to Appendix) numerically based on a six-parameter 

scattering model and capable for different types of waveforms. 

2. The Radar Model (Equation 3.23) 

Using the same form as the scattering model, a similar set of code is built (refer to 
Appendix) numerically based on arbitrary position p = and arbitrary velocity 

u = . The objective will be to determine the values of p and u that maximize 

the cross correlation between the arbitrary model and the measured data. 

3. The Image Model (Equation 3.24) 

The image model forms the main routine of the code (refer to Appendix). It calls 
the parametric radar model as a sub-routine and cross correlated it to the scattering 
model. 
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As part of the development effort to build the full scale imaging algorithm, this 
thesis will focus on imaging a single static point scatterer. The result of the image is also 
used for comparison with the point spread function image. Equation 3.24 does not restrict 
the choice of transmitted waveform nor the transmitter/receiver configuration. The 
objective here is to study the behavior of the new imaging scheme with a single pulse 
waveform and under different transmitter/receiver configuration. 

Figure 8 shows the transmitter/receiver geometries for the simulation performed. 
The receivers are denoted by “A” while the transmitter is denoted by “ + ”. Figure 8(a) 
and (b) are used to show the straight line receiver configuration while Figure 8(c) shows 
the circular receiver configuration. 


8 Receivers 1 Transmitter In Linear(Horizontal) Configuration 
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8 Receivers 1 Transmitter In Linear(Horizontai) Configuration 
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(c) 

Figure 8. Transmitter/Receiver configurations 

A. IMAGING A STATIC POINT TARGET FOR SINGLE RECTANGULAR 
PULSE 

The transmitted waveform used is a single rectangular pulse defined in Equation 
5.6. This waveform is substituted into Equation 3.24 and implemented in Matlab. 
Equation 3.24 is capable of determining a three-dimensional distribution of position in 
space and the corresponding three-dimensional velocities (i.e. a six-dimensional phase 
space image). However, to study the effects of transmitter/receiver geometry and 
translation in position on the imaging scheme, the use of two-dimensional plots of the 
target image in position space will suffice. 

I. Localization of Point-Scatterer in Position Space 

The point target or scatterer is located at position() = (0,0). Two 

dimensional plots of the target image in position space for 4, 9, 36 and 180 receivers are 
given in Eigures 9(a)-(d) respectively. It is shown that the point scatterers are localized in 
position space about the scatterer position (0,0). However, under the different geometries. 
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y-component(km) y-component(km) 

(position) (position) 


the point scatterer is represented differently; this is shown in Figure 10(a)-(d) and 11(a)- 
(d); it can be observed that the geometry of the transmitters/receivers plays an important 
role in shaping the target image. 
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Figure 9. 2-D Localization of point-scatterer in position space (for rectangular pulse). 
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Image of Point Scatterer in Position Space (2-D) 
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Image of Point Scatterer in Position Space (2-D) 
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Figure 10. Transmitter/Receiver in linear (horizontal) configurations 
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Image of Point Scatterer in Position Space (2-D) 
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Figure 11. Transmitter/Receiver in linear (vertical) configurations 


2. Comparison with PSF Image 


When compared with a PSF image, both the images show that the point scatterer 
localized in position space; however it is observed that the PSF gives a smoother spread 
compared to the imaging code. This could be due to the fact that the PSF code is built 
based on an analytical approach where as the imaging code is built based on the 
numerical approach, and as a result there will be some round off errors. 
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Figure 12. Comparison with PSF (target at (v^v^) = (0,0)) 
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B. CONCLUSIONS 


Image artifacts result when motion in the imaging scene is not addressed. In 
recent years, a number of attempts to develop imaging techniques that can handle moving 
objects have been proposed. However, all of these techniques rely on making the 
approximation that a target in motion is assumed to be momentarily stationary while it is 
being interrogated by a radar pulse. This thesis is a study of the new imaging approach 
developed by Cheney and Borden [2] that can accommodate target motion during the 
imaging process. 

In this thesis, the subject of radar imaging from scattered waves is explored and 
applied to the new imaging approach. The simulation result using Matlab obtained 
showed that the new imaging scheme is well behaved; specifically, it localizes the target 
in space and is translation invariant. In addition, it is apparent that the geometry of the 
transmitters and receivers disposition affects the behavior of the imaging system. 

This work represents the early stages of development of a full imaging algorithm 
applicable to moving targets in a multistatic environment. Follow-on work will include 
further simulation efforts. This thesis used single pulse radar waveforms and a single 
transmitter, but in reality, fielded radar systems typically utilize pulse compression 
techniques and send out long, coded pulse trains. Neither need the multistatic radar 
system be restricted to one transmitter. In this case, the imaging scheme has to be 
adjusted to incorporate more realistic radar waveforms and optimize the 
transmitter/receiver geometries for implementing the imaging algorithm. These will be 
used to support the development of the eventual imaging algorithm. Finally, real world 
target data can then be applied to the developed imaging algorithm to assess its 
performance. 
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APPENDIX: MATLAB CODES 


IMAGING CODE 

00000000000000000000000000000000000000000000000000000000000 
%% This Program evaluates the Image of a single point 
target 

% The program uses the approach, techniques on imaging from 
Professor Borden, Brett 

% The Waveform adopted is a 1 x 11 Rect Signal. 

% The Program load the sample_data.mat from the Make_Data.m 
and correlate it to radar model in predict.m subroutine. 
00000000000000000000000000000000000000000000000000000000000 
clear 

load sample_data.mat Rx Ry Tx Ty w k_y T delta_T D N 

N=length(Rx); 

NN=length(Tx); 

J=l; %Set Jacobian to 1 

% Define the area of the Scene (10km by 10km) 

p_x_start=-2 0 ; 

p_x_end=2 0; 

p_y_start=-2 0 ; 

p_y_end=2 0 ; 


res_p=0.5; %Nyquist Sampling rate. 
p_x=p_x_start:res_p:p_x_end; 

P_y=P_y_start:res_p:p_y_end; 
nx=length(p_x); 
ny=length(p_y); 

I=zeros(nx,ny) ; 

for py=l:ny 
for px=l:nx 

P=10''3* [p_x_start+ ( (py-1) *res_p) ;p_y_end- ( (px-1) *res_p) ] ; 
expect=predict(P,Rx,Ry,Tx,Ty,N,NN,w,k_y,T,delta_T); 
for n=l:N 

I(px,py)=I(px,py)+D(n, :)*expect (n, :) '*J; 

end 

end 

end 
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%% Output Plot 
colormap (gray) 

imagesc(p_x,p_y,flipud(abs(I))) 

%colorbar 

xlabel({'\foutsize{12} x-component (km) ' , '\fontsize{12} 

(position)'}) 

ylabel ({'\fontsize{12 } y-component (km) ' , '\fontsize{12} 

(position)'}) 

title('\bf \fontsize{14} \fontname{tahoma} Image of Point 

Scatterer in Position Space (2-D)') 

text(Tx(1)/10^3,Ty(1)/10^3,'\bf \fontsize{14} +') 

for i=l:N 

text(Rx(i)/10^3,Ry(i)/10^3,'\bf \fontsize{9} \Delta') 
end 


SAMPLE DATA MODEL CODE 


00000000000000000000000000000000000000000000000000000000000 
%% This Program Create sample data set 

% The program uses the approach, techniques on imaging from 
Professor Borden, Brett 

00000000000000000000000000000000000000000000000000000000000 


%circular receiver configuration 
r= 10; 

a= 0 :1.5 : 3 6 0; 
b=a*pi/180; 

N=241; 

Rx=10^3*r*cos(b); 

Ry=10^3*r*sin(b); 

% transmitter locations 
M=l; % number of transmitters 
Tx=10^3*[-7]; % x-coords 
Ty=10"^3* [ 0 ] ; % y-coords 

% transmitter locations monostatic case 
%M=N; % number of transmitters 
%Tx=Rx; % x-coords 
%Ty=Ry; % y-coords 

% point targets 
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NN=1; % one target 
XX=10^3*[0]; % x-coords 
YY=10^3*[0]; % y-coords 

% waveform employed (all transmitters use same waveform) 
w=zeros(1,11); 
w(4:8)=l; % step 

%W=fftshift(fft(w)); % sine function 


%%%%%%%%%% make appropriate data 
% Define Variables 


00000000000000 


c=3*10''8; % Speed of Light 

k_y=2*pi*(10^10)/c; % wave vector @f=10 GHz 

%k_y=l ; 

T=2e-3; % Total Data length in sec 


delta_T=T/100; % Data resolution 

D=zeros(N,1000); % large enough to absorb overlap of W and 

the target scene size 


for n=l:N % number of receivers (currently, set up as 
monostatic) 

for nn=l:NN % number of targets 

Rl=sqrt((Tx(1)-XX(nn))^2+(Ty(1)-YY(nn))^2); 

R2 = sqrt( (Rx(n)-XX(nn))^2+(Ry(n)-YY(nn))^2) ; 
tau=(R1+R2)/c; 
m=ceil(tau/delta_T); 

D(n,m+1:m+11)=D(n,m+1:m+11)+w; % superposition principle 

end 

end 


save sample_data.mat Rx Ry Tx Ty w k_y T delta_T D N 
imagesc(abs(D)) 

% read with: load " sample_data. mat" Rx Ry Tx Ty W k_y T 
delta_T D 


PREDICTED DATA MODEL CODE 

00000000000000000000000000000000000000000000000000000000000 
%% This Program creates the subrountine for the radar 
model. 

00000000000000000000000000000000000000000000000000000000000 


o o 
o o 

function d=predict(p,Rx,Ry,Tx,Ty,N,NN,w,k_y,T,delta_T) 
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c=3*10-"8; 
d=zeros(N,1000); 

for n=l:N % number of receivers (currently, set up as 
monostatic) 

Rl = sqrt ( (Tx(l)-p(l) ) . ^2+(Ty (1)-p (2) ) .^2); 

R2 = sqrt ( (Rx (n) -p (1) ) . "^2 + (Ry (n) -p (2) ) . "^2) ; 

tau=mod((R1+R2)/c, 10 0 0); 

m=ceil(tau/delta_T); 

d(n,m+1:m+11)=d(n,m+1:m+11)+w; 

end 


POINT SPREAD FUNCTION CODE 


'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 

% This program evaluates the Point Spread Function (PSF) 

% The program uses the approach, techniques on imaging from 
Professor Borden, % Brett 

% It is a six-parameter scattering model based on arbitrary 
position 

% p=(p_x,p_y,p_z) & velocity u=(u_x,u_y,u_z) 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 

clear all 

%% Variables to be 
x= (10^3)* [0;0;0] ; 

V= [ 0;0;0] ; 

yl= [2;0;0] ; 
zl= [-5;5;0] ; 
z2= [15;15;0] ; 
z3=[-2;-15;0] ; 
z4= [-10;10;0] ; 

Z5=[0;14.1;0] ; 
z6= [10;10;0] ; 

Z7=[14.1;0;0] ; 
z8= [10;-10;0] ; 

U= [ 0;0;0] ; 

pulse=10^(-4); 

%% Parameters Defined 

Yl=(l/sqrt(yl(l)^2+yl(2)^2+yl(3)^2)).*yl; 

Zl=(1/sqrt(zl(1)^2+zl(2)^2+zl(3)^2)).*zl; 

Z2=(1/sqrt(z2(l)^2+z2(2)^2+z2(3)^2)).*z2; 

Z3=(1/sqrt(z3(1)^2+z3(2)^2+z3(3)^2)).*z3; 


specified 

% actual target 
% actual target velocity 
% transmitter_l 
% receiver_l 
% receiver_2 
% receiver 3 
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Z4=(1/sqrt(z4(1)^2+z4 
Z5=(1/sqrt(z5(1)^2+z5 
Z6=(1/sqrt(z6(1)^2+z6 
Z7=(1/sqrt(z7(1)^2+z7 
Z8=(1/sqrt(z8(1)^2+z8 
J=l; 

c=3*10^8; 

k_y=2*pi*(10^10)/c; 

%% PSF Range Defined 
p_x_start=-20; 
p_x_end=2 0; 
p_y_start=-20; 
p_y_end=2 0; 


res p=0. 

1; 



o, 

"o 

resolution of 

P 


p x=p X 
values 

start:res p:p 

X 

end; 

o, 

"o 

define 

vector 

of 

x-axis 

p_y=p_y_ 

start:res p:p 

_y_ 

end; 

o, 

"o 

define 

vector 

of 

y-axis 


values 

%% PSF Velocity Defined 
u_x_start=-200; 
u_x_end=2 00; 
u_y_start=-150; 
u_y_end=150; 
res_u=l; 

u_x=u_x_start:res_u:u_x_end; 
u_y=u_y_start:res_u:u_y_end; 


(2)^2 + z4 (3)^2) ) .*z4; 

(2)^2 + z5 (3)^2) ) .*z5; 

(2)^2 + z6 (3)^2) ) .*z6; 

(2)^2 + z7 (3)^2) ) .*z7; 

(2)^2 + z8 (3)^2) ) .*z8; 

% Jacobian 

% wave vector @f=10 GHz 


'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 

Localization in Position-Space 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 

nx=length(p_x); 
ny=length(p_y); 

K=zeros(ny,nx); 

Kl=zeros(ny,nx); 

K2=zeros(ny,nx); 

K3=zeros(nx,ny); 

K4=zeros(nx,ny); 

K5=zeros(nx,ny); 

K6=zeros(nx,ny); 

K7=zeros(nx,ny); 

K8=zeros(nx,ny); 

for m=l:ny 

for n=l:nx 

p=(10^3)*[p_x_start+((n-1)*res_p);p_y_end-((m- 
1)*res_p);0]; 


% for receiver_l 

argl_l=((Yl+Zl).')*((u*((Z1.')*p))-(v*((Z1.')*x))); 


% length of x = columns in K 
% length of y = rows in K 
% initialize K matrix 
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phasel_l=exp(-i*k_y*argl_l/c); 
arg2_l=((Yl+Zl).')*((p-x)+(0.5*(u- 
v)*(((Yl+Zl).')*(p+x)))); 

phase2_l=exp(i*k_y*arg2_l/c); 
tau_l=l/c*(((Yl+Zl).')*(p-x)); 
neu_l=k_y*(((Yl+Zl).')*(u-v)); 

%rectangular pulse 
amb_l=pi*(pulse-abs(tau_l))*exp(- 
i*0.5*neu_l*pulse)*sinc((0.5/pi)*(-neu_l)*(pulse- 
abs (tau_l))); 

K1(m,n)=phasel_l*phase2_l*amb_l*J; 

% for receiver_2 

argl_2=((Y1+Z2).')*((u*((Z2.')*p))-(v*((Z2.')*x))) 
phasel_2=exp(-i*k_y*argl_2/c); 
arg2_2=((Yl+Zl).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z2).')*(p+x)))); 

phase2_2=exp(i*k_y*arg2_2/c); 
tau_2=l/c*(((Y1+Z2).')*(p-x)); 
neu_2=k_y*(((Y1+Z2).')*(u-v)); 
amb_2=pi*(pulse-abs(tau_2))*exp(- 
i*0.5*neu_2*pulse)*sinc((0.5/pi)*(-neu_2)*(pulse- 
abs (tau_2))); 

K2(m,n)=phasel_2*phase2_2*amb_2*J; 

% for receiver_3 

argl_3=((Y1+Z3).')*((u*((Z3.')*p))-(v*((Z3.')*x))) 
phasel_3=exp(-i*k_y*argl_3/c); 
arg2_3=((Y1+Z3).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z3).')*(p+x)))); 

phase2_3=exp(i*k_y*arg2_3/c); 
tau_3=l/c*(((Y1+Z3).')*(p-x)); 
neu_3=k_y*(((Y1+Z3).')*(u-v)); 
amb_3=pi*(pulse-abs(tau_3))*exp(- 
i*0.5*neu_3*pulse)*sinc((0.5/pi)*(-neu_3)*(pulse- 
abs (tau_3))); 

K3(m,n)=phasel_3*phase2_3*amb_3*J; 

% for receiver_4 

argl_4=((Y1+Z4).')*((u*((Z4.')*p))-(v*((Z4.')*x))) 
phasel_4=exp(-i*k_y*argl_4/c); 
arg2_4=((Y1+Z4).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z4).')*(p+x)))); 

phase2_4=exp(i*k_y*arg2_4/c); 
tau_4=l/c*(((Y1+Z4).')*(p-x)); 
neu_4=k_y*(((Y1+Z4).')*(u-v)); 
amb_4=pi*(pulse-abs(tau_4))*exp(- 
i*0.5*neu_4*pulse)*sinc((0.5/pi)*(-neu_4)*(pulse- 
abs (tau_4))); 

K4(m,n)=phasel_4*phase2_4*amb_4*J; 

% for receiver_5 

argl_5=((Y1+Z5).')*((u*((Z5.')*p))-(v*((Z5.')*x))) 
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phasel_5=exp(-i*k_y*argl_5/c); 
arg2_5=((Y1+Z5).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z5).')*(p+x)))); 

phase2_5=exp(i*k_y*arg2_5/c); 
tau_5=l/c*(((Y1+Z5).')*(p-x)); 
neu_5=k_y*(((Y1+Z5).')*(u-v)); 
amb_5=pi*(pulse-abs(tau_5))*exp(- 
i*0.5*neu_5*pulse)*sinc((0.5/pi)*(-neu_5)*(pulse- 
abs (tau_5))); 

K5(m,n)=phasel_5*phase2_5*amb_5*J; 

% for receiver_6 

argl_6=((Y1+Z6).')*((u*((Z6.')*p))-(v*((Z6.')*x))); 
phasel_6=exp(-i*k_y*argl_6/c); 
arg2_6=((Y1+Z6).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z6).')*(p+x)))); 

phase2_6=exp(i*k_y*arg2_6/c); 
tau_6=l/c*(((Y1+Z6).')*(p-x)); 
neu_6=k_y*(((Y1+Z6).')*(u-v)); 
amb_6=pi*(pulse-abs(tau_6))*exp(- 
i*0.5*neu_6*pulse)*sinc((0.5/pi)*(-neu_6)*(pulse- 
abs (tau_6))); 

K6(m,n)=phasel_6*phase2_6*amb_6*J; 

% for receiver_7 

argl_7=((Y1+Z7).')*((u*((Z7.')*p))-(v*((Z7.')*x))); 
phasel_7=exp(-i*k_y*argl_7/c); 
arg2_7=((Y1+Z7).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z7).')*(p+x)))); 

phase2_7=exp(i*k_y*arg2_7/c); 
tau_7=l/c*(((Y1+Z7).')*(p-x)); 
neu_7=k_y*(((Y1+Z7).')*(u-v)); 
amb_7=pi*(pulse-abs(tau_7))*exp(- 
i*0.5*neu_7*pulse)*sinc((0.5/pi)*(-neu_7)*(pulse- 
abs (tau_7))); 

K7(m,n)=phasel_7*phase2_7*amb_7*J; 

% for receiver_8 

argl_8=((Y1+Z8).')*((u*((Z8.')*p))-(v*((Z8.')*x))); 
phasel_8=exp(-i*k_y*argl_8/c); 
arg2_8=((Y1+Z8).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z8).')*(p+x)))); 

phase2_8=exp(i*k_y*arg2_8/c); 
tau_8=l/c*(((Y1+Z8).')*(p-x)); 
neu_8=k_y*(((Y1+Z8).')*(u-v)); 
amb_8=pi*(pulse-abs(tau_8))*exp(- 
i*0.5*neu_8*pulse)*sinc((0.5/pi)*(-neu_8)*(pulse- 
abs (tau_8))); 

K8(m,n)=phasel_8*phase2_8*amb_8*J; 

K(m,n)=(abs(K1(m,n)+K2(m,n)+K3(m,n)+K4(m,n)+K5(m,n)+K6(m,n)+ 
K7(m,n)+K8(m,n))).^2; 
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end 


end 

9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 

'o'o'o'o'o'o'o'o'o'o o'o'o'o'o'o'o 

Localization in Velocity-Space 

9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 

nx=length(u_x); % length of x = columns in K 

ny=length(u_y); % length of y = rows in K 

K=zeros(ny,nx); 

Kl=zeros(ny,nx); 

K2=zeros(ny,nx); 

K3=zeros(ny,nx); 

K4=zeros(ny,nx); 

K5=zeros(ny,nx); 

K6=zeros(ny,nx); 

K7=zeros(ny,nx); 

K8=zeros(ny,nx); 

for m=l:ny 

for n=l:nx 

u=[u_x_start+((n-1)*res_u);u_y_end-((m-1)*res_u);0]; 
% for receiver_l 

argl_l=((Yl+Zl).')*((u*((Z1.')*p))-(v*((Z1.')*x))); 
phasel_l=exp(-i*k_y*argl_l/c); 
arg2_l=((Yl+Zl).')*((p-x)+(0.5*(u- 
v)*(((Yl+Zl).')*(p+x)))); 

phase2_l=exp(i*k_y*arg2_l/c); 
tau_l=l/c*(((Yl+Zl).')*(p-x)); 
neu_l=k_y*(((Yl+Zl).')*(u-v)); 
amb_l=pi*(pulse-abs(tau_l))*exp(- 
i*0.5*neu_l*pulse)*sinc((0.5/pi)*(-neu_l)*(pulse- 
abs (tau_l))); 

K1(m,n)=phasel_l*phase2_l*amb_l*J; 

% for receiver_2 

argl_2=((Y1+Z2).')*((u*((Z2.')*p))-(v*((Z2.')*x))); 
phasel_2=exp(-i*k_y*argl_2/c); 
arg2_2=((Yl+Zl).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z2).')*(p+x)))); 

phase2_2=exp(i*k_y*arg2_2/c); 
tau_2=l/c*(((Y1+Z2).')*(p-x)); 
neu_2=k_y*(((Y1+Z2).')*(u-v)); 
amb_2=pi*(pulse-abs(tau_2))*exp(- 
i*0.5*neu_2*pulse)*sinc((0.5/pi)*(-neu_2)*(pulse- 
abs (tau_2))); 

K2(m,n)=phasel_2*phase2_2*amb_2*J; 

% for receiver_3 

argl_3=((Y1+Z3).')*((u*((Z3.')*p))-(v*((Z3.')*x))); 
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phasel_3=exp(-i*k_y*argl_3/c); 
arg2_3=((Y1+Z3).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z3).')*(p+x)))); 

phase2_3=exp(i*k_y*arg2_3/c); 
tau_3=l/c*(((Y1+Z3).')*(p-x)); 
neu_3=k_y*(((Y1+Z3).')*(u-v)); 
amb_3=pi*(pulse-abs(tau_3))*exp(- 
i*0.5*neu_3*pulse)*sinc((0.5/pi)*(-neu_3)*(pulse- 
abs (tau_3))); 

K3(m,n)=phasel_3*phase2_3*amb_3*J; 

% for receiver_4 

argl_4=((Y1+Z4).')*((u*((Z4.')*p))-(v*((Z4.')*x))) 
phasel_4=exp(-i*k_y*argl_4/c); 
arg2_4=((Y1+Z4).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z4).')*(p+x)))); 

phase2_4=exp(i*k_y*arg2_4/c); 
tau_4=l/c*(((Y1+Z4).')*(p-x)); 
neu_4=k_y*(((Y1+Z4).')*(u-v)); 
amb_4=pi*(pulse-abs(tau_4))*exp(- 
i*0.5*neu_4*pulse)*sinc((0.5/pi)*(-neu_4)*(pulse- 
abs (tau_4))); 

K4(m,n)=phasel_4*phase2_4*amb_4*J; 

% for receiver_5 

argl_5=((Y1+Z5).')*((u*((Z5.')*p))-(v*((Z5.')*x))) 
phasel_5=exp(-i*k_y*argl_5/c); 
arg2_5=((Y1+Z5).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z5).')*(p+x)))); 

phase2_5=exp(i*k_y*arg2_5/c); 
tau_5=l/c*(((Y1+Z5).')*(p-x)); 
neu_5=k_y*(((Y1+Z5).')*(u-v)); 
amb_5=pi*(pulse-abs(tau_5))*exp(- 
i*0.5*neu_5*pulse)*sinc((0.5/pi)*(-neu_5)*(pulse- 
abs (tau_5))); 

K5(m,n)=phasel_5*phase2_5*amb_5*J; 

% for receiver_6 

argl_6=((Y1+Z6).')*((u*((Z6.')*p))-(v*((Z6.')*x))) 
phasel_6=exp(-i*k_y*argl_6/c); 
arg2_6=((Y1+Z6).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z6).')*(p+x)))); 

phase2_6=exp(i*k_y*arg2_6/c); 
tau_6=l/c*(((Y1+Z6).')*(p-x)); 
neu_6=k_y*(((Y1+Z6).')*(u-v)); 
amb_6=pi*(pulse-abs(tau_6))*exp(- 
i*0.5*neu_6*pulse)*sinc((0.5/pi)*(-neu_6)*(pulse- 
abs (tau_6))); 

K6(m,n)=phasel_6*phase2_6*amb_6*J; 

% for receiver_7 

argl_7=((Y1+Z7).')*((u*((Z7.')*p))-(v*((Z7.')*x))) 
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phasel_7=exp(-i*k_y*argl_7/c); 
arg2_7=((Y1+Z7).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z7).')*(p+x)))); 

phase2_7=exp(i*k_y*arg2_7/c); 
tau_7=l/c*(((Y1+Z7).')*(p-x)); 
neu_7=k_y*(((Y1+Z7).')*(u-v)); 
amb_7=pi*(pulse-abs(tau_7))*exp(- 
i*0.5*neu_7*pulse)*sinc((0.5/pi)*(-neu_7)*(pulse- 
abs (tau_7))); 

K7(m,n)=phasel_7*phase2_7*amb_7*J; 

% for receiver_8 

argl_8=((Y1+Z8).')*((u*((Z8.')*p))-(v*((Z8.')*x))); 
phasel_8=exp(-i*k_y*argl_8/c); 
arg2_8=((Y1+Z8).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z8).')*(p+x)))); 

phase2_8=exp(i*k_y*arg2_8/c); 
tau_8=l/c*(((Y1+Z8).')*(p-x)); 
neu_8=k_y*(((Y1+Z8).')*(u-v)); 
amb_8=pi*(pulse-abs(tau_8))*exp(- 
i*0.5*neu_8*pulse)*sinc((0.5/pi)*(-neu_8)*(pulse- 
abs (tau_8))); 

K8(m,n)=phasel_8*phase2_8*amb_8*J; 

K(m,n)=(abs(K1(m,n)+K2(m,n)+K3(m,n)+K4(m,n)+K5(m,n)+K6(m,n)+ 
K7(m,n)+K8(m,n))).^2; 

end 

end 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 

3-D Slice Plot for Position-Space 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 


p_x_start=-20; 
p_x_end=2 0; 
p_y_start=-20; 
p_y_end=2 0; 
p_z_start=-20; 
p_z_end=2 0; 
res_p=0.5; 

p_x=p_x_start:res_p:p_x_end; 
values 

p_y=p_y_start:res_p:p_y_end; 
values 

p_z=p_z_start:res_p:p_z_end; 

nx=length(p_x); 
ny=length(p_y); 
nz=length(p_z); 

K=zeros(ny,nx,nz); 

Kl=zeros(ny,nx,nz); 


% define vector of x-axis 
% define vector of y-axis 


% length of x = columns in K 

% length of y = rows in K 

% length of z = pages in K 

% initialize K matrix 
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K2=zeros(ny,nx,nz); 


for a=l:nz 

for m=l:ny 

for n=l:nx 

p=(10^3)*[p_x_start+((n-1)*res_p);p_y_end-((m- 
1)*res_p);p_z_start+((a-1)*res_p)]; 

% for receiver_l 

argl_l=((Yl+Zl).')*((u*((Z1.')*p))- 
(v*((Zl.')*x))); 

phasel_l=exp(-i*k_y*argl_l/c); 
arg2_l=((Yl+Zl).')*((p-x)+(0.5*(u- 
v)*(((Yl+Zl).')*(p+x)))); 

phase2_l=exp(i*k_y*arg2_l/c); 
tau_l=l/c*(((Yl+Zl).')*(p-x)); 
neu_l=k_y*(((Yl+Zl).')*(u-v)); 
amb_l=pi*(pulse-abs(tau_l))*exp(- 
i*0.5*neu_l*pulse)*sinc((0.5/pi)*(-neu_l)*(pulse- 
abs (tau_l))); 

K1(m,n,a)=phasel_l*phase2_l*amb_l*J; 

% for receiver_2 

argl_2=((Y1+Z2).')*((u*((Z2.')*p))- 
(v*((Z2.')*x))); 

phasel_2=exp(-i*k_y*argl_2/c); 
arg2_2=((Yl+Zl).')*((p-x)+(0.5*(u- 
v)*(((Y1+Z2).')*(p+x)))); 

phase2_2=exp(i*k_y*arg2_2/c); 
tau_2=l/c*(((Y1+Z2).')*(p-x)); 
neu_2=k_y*(((Y1+Z2).')*(u-v)); 
amb_2=pi*(pulse-abs(tau_2))*exp(- 
i*0.5*neu_2*pulse)*sinc((0.5/pi)*(-neu_2)*(pulse- 
abs (tau_2))); 

K2(m,n,a)=phasel_2*phase2_2*amb_2*J; 

K(m,n,a)=(abs(K1(m,n,a)+K2(m,n,a))).^2; 


end 

end 

end 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 

Ambiguity Implementation for Chirp Signal 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 


amb_l=pi*(pulse-abs(tau_l))*exp(i*0.5*(- 
neu_l+(gamma*tau_l))*pulse)*sinc((0.5/pi)*(- 
neu_l+(gamma*tau_l))*(pulse-abs(tau_l))); 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 


71 


Output Plot 

9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 9 ' 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 


clear all 

image(p_x,p_y,64*K/max(max(K))) 
hold on 

plot(2,0,'+k',-5,5,'^k',15,15,'^k',-2,-15,'^k',15,- 
15,'^k','linewidth',2) 
hold off 

xlabel({'\fontsize{11} x-component','\fontsize{11} 
(position)'}) 

ylabel({'\fontsize{11} y-component','\fontsize{11} 
(position)'}) 

title({'\fontsize{12} Transmitter/Receiver Geometry 1'}) 

% Slice Plot 
norm=max(max(max(K))); 

[TJX, UY, UZ] =meshgrid (u_x, u_y, u_z) ; 
slice(UX,UY,UZ,K/norm,[0],[0],[0]) 
colorbar 
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